On The Existence of Roton Excitations in Bose Einstein Condensates: Signature of 

Proximity to a Mott Insulating Phase 
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Within the last decade, artificially engineered Bose Einstein Condensation has been achieved in 
atomic systems. Bose Einstein Condensates are superfluids just like bosonic Helium is and all 
interacting bosonic fluids are expected to be at low enough temperatures. One difference between 
the two systems is that superfluid Helium exhibits roton excitations while Bose Einstein Condensates 
have never been observed to have such excitations. The reason for the roton minimum in Helium is 
its proximity to a solid phase. The roton minimum is a consequence of enhanced density fluctuations 
at the reciprocal lattice vector of the stillborn solid. Bose Einstein Condensates in atomic traps 
are not near a solid phase and therefore do not exhibit roton minimum. We conclude that if Bose 
Einstein Condensates in an optical lattice are tuned near a transition to a Mott insulating phase, a 
roton minimum will develop at a reciprocal lattice vector of the lattice. Equivalently, a peak in the 
structure factor will appear at such a wavevector. The smallness of the roton gap or the largeness 
of the structure factor peak are experimental signatures of the proximity to the Mott transition. 



In the present work we focus attention in the possible 
existence or nonexistence of roton excitations in artifi- 
cially engineered Bose Einstein Condensates (BECs)0, 
0, HI . BECs are superfluid as they posses a finite sound 
speed 0,01 El an d reduced long wavelength scattering 
The other important bosonic superfluid is He4 0,] 
Since BECs and He4 are the same universal phase of 
matter, they share a large commonality. For example, 
they are both superfluids whose elementary long wave- 
length excitation spectrum is phononic, their ground 
states spontaneously break gauge invariance[T3, that is 
their ground states exhibit macroscopic quantum coher- 
ence, and they both posses vortex excitations^!). This is 
an example of universality of stable fixed points of mat- 
ter. In these cases the superfluid ground state is such 
a fixed point. The similarities will become ever more 
apparent as the low energy properties BECs are studied 
more carefully. 

On the other hand, at shorter wavelengths Helium pos- 
sesses a roton minimum|l i 11^ ] in its excitation spec- 
trum leading to low energy excitations occurring at a 
specific wavevector, which are gapped and are different 
than long wavelength sound. These excitations can make 
important contributions to the dynamics and thermody- 
namics of the system. In BECs no such excitation has 
been observed and the lack of a peak in the structure 
factor leads to the conclusion that they do not posses 
roton excitations. We review some of the properties of 
rotons in Helium and determine under what conditions 
will they occur in BECs. 

All quantum many particle systems are described by a 
Hamiltonian 5] 



H = 



v ■ pv 



d 3 r + U[p] 



(1) 



where mv ■ pv/2 is the kinetic energy operator, U[p] is 
the potential energy operator, which in general can be a 



functional of the density operator. The density operator 
in first quantized notation is 



(2) 



with fj being the position of the ith particle. The velocity 
operator is in first quantized notation 
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with commutation relations 



[p(r ') ,v(r)] 



i-V5(f' 
m 



f) 



(3) 



(4) 



The ground state of the Hamiltonian could be a quan- 
tum fluid or solid depending on the interaction. We will 
suppose it to be a fluid as we have BECs and Helium in 
mind. 

We study the excitation spectrum of the bosonic fluid 
on quite general grounds following LandauQ and the 
Russian school 01 quite closely. We suppose the den- 
sity to have a well defined average which we take to be 
uniform for simplicity. This will not change the nature 
of the considerations, although in real life the density 
can be modulated by the lattice as for BECs in optical 
lattices. The average velocity is zero as we consider a sys- 
tem at rest. We Fourier expand the density and velocity 
operators about their averages: 



1 



P(r) = Po + ^ 



ik-r 



(5) 



(6) 



where we are using the lattice normalization, the number 
of sites N instead of the volume V. Quantum mechani- 
cally, the velocity is proportional to the gradient of the 
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local phase variable of the particle, v(f) = H/mVd, or in 
Fourier components 



ihk 

v k = e k 

m 



(7) 



The density velocity commutation relation (J2J imply the 
well known density phase commutation relation 



■s>] = % kw 



(8) 



i.e. hQ_k is the momentum conjugate to the density fluc- 
tuation ,0£. The Hamiltonian thus becomes 



rn 



n 2 k 2 



±k ~ p k 2 ' JV; < ' ' ' < (10) 

The density energy excitation spectrum of the quantum 
liquid is 



Er = hujAn + 1/2) 



2 K k 



Pok 2 + 

k impo 



(11) 



(12) 



The ground state energy is given by U(p ) + ^2^huj^/2. 
We note that the Virial theorem implies that 



2 k 



2 1,2 



^k 



2mS(k) 



(13) 



where the structure factor is defined by 

(\Pk\ 2 ) 



K 

(9) 

where U? is the Fourier transform of the interaction. The 
interaction is completely general. The only requirements 
are that Uz~ is a constant, so that we can have a super- 
fluid, and that there is some short range repulsion analo- 
gous to a Hubbard term that can be turned up in order to 
stabilize a Mott insulating phase, (h 2 k 2 ) / (4mpo)\pj;\ 2 is 
the elastic energy associated with changing the density. 
This last term comes from acting the Hamiltonian on 
the wavefunction and making sure that one keeps careful 
track of both density and phase degrees of freedom [T^ 
in the full quantum mechanical kinetic energy term. We 
have performed studies very near to these in the context 
of fermionic systems) 15$ which we follow very closely in 
the present note. Charged boson systems were studied 
previously^] with the Quantum Hall effect in mind. 

The quantum liquid is thus a collection of harmonic os- 
cillators in momentum space for the density fluctuations. 
The mass and spring constant of the oscillators are 



S(k) 



Np 



(14) 



which is, of course, the Fourier transform of the density- 
density correlation function. For the quantum liquid we 
thus have 



S(k) 



hk 



l 



2m 



jt I h 2 k 2 
u k' Ampa 



(15) 



The ground state energy will have the density oscillators 
unexcited. Hence, if the ground state wavefunction of 
the bosonic system without density correlations is \ij)Gs), 
the ground state wavefunction with density correlations 



exp 



poh 2 k 2 



h 2 k 2 

Ampo 



i>A 2 )VPgs) 
(16) 

The factor in front appears because it is the ground state 
wavefunction of harmonic oscillators and the minimum 
energy condition is obviously to have the oscillators shak- 
ing as little as possible. We see that the factor suppresses 
density fluctuations. The wavefunction is essentially ex- 
act to quadratic order in the density fluctuations. 

Since the Hamiltonian and ground state wavefunction 
do not contain any periodic lattice effects, they aptly de- 
scribe systems with no external lattice such as liquid He- 
lium, to which we turn our discussion now. The ground 
state wavefunction (|l(j|l does not contain any insulating 
crystalline order as (p%) = for all k. When the sys- 
tem spontaneously solidifies, it will exhibit Bragg peaks 
at the reciprocal lattice vectors, Q, corresponding to the 
insulating crystal. These peaks mean that (Pq) 0. 
We inquire into the behavior when the boson fluid ap- 
proaches a solidification transition. 

Proximity to a solid phase with reciprocal lattice vec- 
tor Q for the incipient order softens the spring constant 
of the density oscillators of the liquid, — > 0, which 
makes to^ — ► 0. This last condition means the mode be- 
comes degenerate with the ground state and there will 
be enhanced density fluctuations, — > oo at the 
wavevector of the stillborn solid, i.e. the system is ap- 
proaching a phase transition. The diverging fluctuations 
follow because the position fluctuations of a harmonic os- 
cillator diverge as the spring constant vanishes due to the 
Heisenberg uncertainty principle. The softening of the 
density mode implies that the interaction Ut has a mini- 
mum at Q which in turn produces a peak in the structure 
factor, S(Q), that diverges as one approaches the point 
when the system solidifies. This peak in the structure 
factor implies will drag the spectrum down according to 
the equation (|13fl creating a roton minimum |l2|. 
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FIG. 1: Structure factor vs. wavevector for a Bose fluid near 
a solidification transition. The density fluctuation peak at the 
wavevector of the incipient order is clearly visible. 



the wells of the optical lattice are made deep, the system 
is far from a solid or Mott insulating phase. In fact, it is 
an extremely good approximation to take the interaction 
potential to be a constant U? = I. For such a case, we 
plot in Figures 3 and 4 the structure factor and density 
excitation spectrum. We see no peak in the structure 
factor or roton minimum in the excitation spectrum as 
observed experimentally 7]. 
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FIG. 2: Frequency dispersion for a Bose fluid near a solidifi- 
cation transition. The roton minimum at the wavevector of 
the incipient order is clearly visible. 



We thus have as a simple phenomenological model that 
captures the physics 



U 



u 



2 1.2 



k 2 



k 



k 2 



■{k-Qf-a 



K 2 k 



-A 



where fco is some microscopic scale that characterizes 

the interaction, nonzero a makes a roton minimum, and 

h 2 o 2 

A — is the roton gap, which will collapse at the 

crystallization transition as the smallness of the gap of 
such a roton is a measure of proximity to crystallization. 
The first term in the potential represents the short range 
repulsion, chosen for simplicity to have a Yukawa form. 
The parameter A can be thought of as the long range 
repulsion. From the form of the interaction we get the 
structure factor and the density excitation spectrum for 
Helium, which we plot in Figures 1 and 2. We emphasize 
that the reason for the roton minimum in the excitation 
spectrum is the proximity of the liquid Helium ground 
state to a solid phase. Helium is an almost solid barely 
melted by quantum fluctuations, i.e. it is a very strongly 
correlated liquid. 

We now turn our attention to BECs in optical lattices. 
These superfluids do not exhibit roton minimum in their 
excitation spectrum. The reason for this is that unless 



FIG. 3: Structure factor vs. wavevector for a Bose fluid far 
from a solidification transition. There is no peak. 




FIG. 4: Frequency dispersion for a Bose fluid far from a solid- 
ification transition. The upper graph is the superfluid disper- 
sion and the lower one is a phononic linear dispersion in order 
to make aparent the linearity of the superfluid dispersion at 
long wavelengths. There is no roton minimum. 

So far we have neglected optical lattice effects because 
even though they make quantitative changes to mea- 
surable parameters of the BEC superfluid, they do not 
change the universal superfluid physics. On the other 
hand, as one approaches a commensurate Mott insulat- 
ing phase, the lattice will become ever more important. 
The presence of the optical lattice implies that there is 
an external potential 
adding to the Hamiltonian (JHJ the term 



Vt that acts on the bosons thus 



H ext = - V kPk 



(17) 



The effect of the extra term in the ground state wavefunc- 
tion is to move the equilibrium position of the density 
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oscillators from (p£) = to 



h?k 2 



(18) 



The fact that there is (p?) ^ does not mean the bo- 
son fluid in an optical lattice has solidified into an insu- 
lating phase. It just means that the boson fluid is being 
modulated by the periodic lattice. It is true that in the 
insulating phase this expected value will be nonzero too. 
If there was not an external lattice, p^ ^ for some fc's 
would imply an insulating phase as the only way this can 
happen is for the system to spontaneously break trans- 
lational invariance. So the question of how to determine 
if one is in the Mott phase is delicate for lattice boson 
systems. 
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Roton Gap 

FIG. 5: Density Fourier component at the wavevector of the 
incipient order vs. roton gap for a Bose fluid near a solidifica- 
tion transition. The divergence as the roton gap goes to zero 
is a signature of the imminent transition. 



a large similarity of intrinsically quantum behavior for 
their ground states and low energy states. One impor- 
tant difference is that Helium possesses roton excitations 
while BECs appear not to have them. The reason for 
this is that BECs in the atomic traps are far from being 
a solid. Helium posses a roton because it is a stillborn 
solid where quantum fluctuations barely prevent solidifi- 
cation. 

On the other hand, BECs in optical lattices 0,0,0 
can be tuned near a solidification or Mott transition by 
making the wells deeper. We thus predict that near a 
transition from a superfluid BEC to a commensurate 
Mott phase in an optical lattice there will be roton ex- 
citations measurable by light scattering occurring at the 
reciprocal lattice vector of the lattice. Their experimen- 
tal signatures will be a large peak in the structure factor 
or a roton minimum in the excitation spectrum. The 
smallness of the roton gap or height of the structure fac- 
tor peak can be used as diagnostics of the proximity to 
the incipient Mott order. The roton excitation adds to 
the list of intrinsically quantum mechanical excitations 
that exists in BECs and will be of importance to the 
dynamics of the systems. 
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A way to know that the bosons in an optical lattice 
have solidified into a commensurate Mott phase is to look 
at the density excitation spectrum or the structure factor. 
When it has solidified, the excitation spectrum will have 
a gap at long wavelengths and the structure factor will 
go like k 2 at long wavelengths. This is not a violation of 
Goldstone's theorem^] as the lattice was not generated 
spontaneously by the boson system but is provided by the 
external lattice potential, i.e. the symmetry is explicitely 
broken. Our interest is what are the signatures of the in- 
cipient Mott insulating phase. As mentioned above, one 
can know the transition is approaching by a softening 
of the excitation spectrum at appropriate wavevectors. 
In this case, the system will behave exactly as Helium, 
exhibiting roton minimum and its structure factor and 
energy excitation spectrum will be as shown in Figures 
1 and 2. Moreover, from the formula for the response of 
the lattice potential we see that the density Fourier 
component at the roton wavevector will grow as the ro- 
ton gap collapses. We plot the growth of such a Fourier 
component vs. size of the roton gap in Figure 5. 

Summarizing, artificially engineered BECs and liquid 
He4 are both quantum bosonic superfluids. They share 
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